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ON THE DEVELOPMENT OF TURBULENT WAKES FROM VORTEX STREETS ‘ 

By Anatol Roshko 



SUMMARY , 

Wake development behind circular cylinders at Reynolds 
numbers Jrom 40 to 10,000 was investigated in a low-speed wind 
tunnel. Standard hot-wire techniques were used to study the 
velocity fluctuations. . 

The Reynolds number range of periodic vortex shedding is 
divided into two distinct subranges. At R=Jfl to 150, called 
the stable range, regular vortex streets are formed and no turbu- 
lent motion is developed. The range R=150 to 300 is a 
transition range to a regime called the irregular range, in which 
turbulent velocity fluctuations accompany the periodic formation 
of vortices. The turbulence is initiated by laminar-turbulent 
transition in the free layers which spring from Jhe separation 
points on the cylinder. This transition first occurs in the 
range R=160 to 300, 

Spectrum and statistical measurements were made to study 
the velocity fluctuations. In the stable range the vortices decay 
by viscous diffusion. In the irregular range the diffusion is 
turbulent and the wake becomes fully turbulent in 40 to 50 
diameters downstream. 

It was found that in the stable range the vortex street has a 
periodic spanwise structure. 

The dependence of shedding freguency on velocity was success- 
fully used to measure flow velocity. 

Measurements in the wake of a ring showed that an annular 
vortex street is developed. 

INTRODUCTION 

It is alw^ays difficult to determine precisely the date and 
author of a discovery or idea. This seems to be the case 
with the periodic phenomena associated with flow about a 
cylinder. Although the effect of wind in producing vibra- 
tions in wires (aeolian tones) had been known for some time, 
the first experimental observations are due to Strouhal (ref. 
1) who showed that the frequency depends on the relative 
air velocity and not the elastic properties of the wires. Soon 
after, Rayleigh (1879, refs. 2 and 3) performed similar 
experiments. His formulation of the Reynolds numbei 
dependence demonstrates his remarkable insight into the 
problem. 

These earliest oliservations were concerned with the rela- 
tions between vibration frequency and wind velocity. The 
periodic nature of the wake was discovered later, although 
Leonardo da Vinci in the fifteenth century had already drawn 
some rather a(*curate sketclu's of the vort(*x formation in the 
flow behind l)luff bodies (ref. 4). However, Leonardo's 


drawings show a symmetric row of vortices in the wake. 
The first modern pictures showing the alternating arrange- 
ment of vortices in the wake were published by Ahlborn in 
1902 (ref. 5); his visualization techniques have been used 
extensively since then. The importance of this phenomenon, 
now known as the Kdrman vortex street, was pointed out by 
Benard (1908, ref. 6). 

In 1911 KArmdn gave his famous theory of the vortex 
street (ref. 7), stimulating a widespread and lasting series of 
investigations of the subject. For the most part these con- 
cerned themselves with experimental comparisons of real 
vortex streets with Kdrmdn^s idealized model, calculations 
on the effects of various disturbances and configurations, and 
so on. It can hardly be said that any fundamental advance 
in the problem has been made since Kdrmdn^s stability 
papers, in which he also clearly outlined the natm*e of the 
phenomenon and the unsolved problems. Outstanding per- 
haps is the problem of the periodic vortex-shedding mecha- 
nism, for which there is yet no suitable theoretical treatment. 

However, the results of the many vortex-street studies, 
especially the experimental ones, are very useful for further 
progress in the problem. Attention should be drawn to the 
work of Fage and his associates (1927, refs. 8 to 10), whose 
experimental investigations were conducted at Reynolds 
numbem well above the ranges examined by most other 
investigators. Their measurements in the wake close behind 
a cylinder provide much useful information about the nature 
of the shedding. More recently Kovasznay (1949, ref. 11) 
has conducted a hot-wire investigation of the stable vortex 
street (low Reynolds numbei-s), to which frequent reference 
will be made. 

Vortex-street patterns which are stable and well defined 
for long distances downstream actually occur in only a small 
range of cylinder Reynolds numbers, from about /? = 40 to 
150, and it is to this range that most of the attention has 
been given. On the other hand, as is well known, periodic- 
vortex shedding also occui*s at higher Reynolds numbei*s, 
up to 10^ or more, but the free vortices which move down- 
stream are quickly obliterated, by turbulent diffusion, and a 
turbulent wake is established. 

The present interest in the vortex street is due to some 
(piestions arising from the study of turbulent flow behind 
cylinders and grids. Such studies are usually made at Reyn- 
olds uumbei*s for which periodic vortex shedding from the 
cylinders or grid rods might occur. However, the measure- 


ments are always taken downstream far enough to insure 
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that the periodic velocity fluctuations are obliterated and 
the flow is completely turbulent. There are several 
important consequences of this limitation. 

First, the energy of the velocity fluctuations is quite low 
compared with the energy near the cylinder, and especially 
low compared with the dissipation represented by the form 
drag. In attaining the developed downstream state there is 
evidently not only a rapid redistribution of energy among 
the spectral components but also a large dissipation. Sec- 
ond, the theories which describe these downstream stages do 
not relate the flow to the initial conditions except very loosely 
in terms of dimensionless parameters, and it is usually 
necessary to determine an origin empirically (e. g., mixing- 
length theory or similarity theories). 

On the other hand, there is evidence that some features 
are permanent, so that they must be determined near the 
beginning of the motion. One such feature is the low-wave- 
number end of the spectrum which (in the theory of homo- 
geneous tiu-bulence) is invariant. 

Another is the random element. It has been pointed out 
by Dryden (refs. 12 and 13) that in the early stages of the 
decay of isotropic turbulence behind grids the bulk of the 
turbulent energy lies in a spectral range which is well ap- 
proximated by the simple function characteristic 

of certain random processes. Liepmann (ref. 14) has sug- 
gested that such a random process may be found in the 
shedding of vortices from the grids. 

In short, there has been no description, other than very 
qualitative,_of the downstream development of wakes which, 
over a wide range of Reynolds number, exhibit a definite 
periodicity at the beginning. The measurements reported 
here were undertaken to help bridge this gap. 

The main results show the downstream development of the 
wake, in terms of energy, spectrum, and statistical properties. 
This development is quite different in two Reynolds number 
ranges, the lower one extending from about 40 to 150 and the 
upper, from 300 to 10< (and probably 10»), with a transition 
range between. The lower range is the region of the classic 
vortex street, stable and regular for a long distance down- 
stream. The fluctuating energy of the flow has a discrete 
spectrum and simply decays downstream without transfer of 
energy to other frequencies. Irregular fluctuations are not 
developed. In the upper range there is still a predominant 
(shedding) frequency in the velocity fluctuations near the 
cylinder, and most of the energy is concentrated at this 
frequency; however, some irregularity is already developed, 
and this corresponds to a continuous spectral distribution of 
some of the energy. Downstream, the discrete energy, at 
the shedding frequency, is quickly dissipated or transferred 
to other frequencies, so that by 50 diameters the wake is 
completely turbulent, and the energy spectrum of the velocity 
fluctuations approaches that of isotropic turbulence 
All other features of the periodic shedding and wake phe- 
nomena may be classified as belonging to one or the other of 
Uie two ranges. This viewpoint allows some systematization 
in the study of wake development. 


In particular, it is felt the possibilities of the vortex street 
are by no means exhausted. A study of the interaction of 
periodic fluctuations with a turbulent fleld seems to be a 
fruitful approach to the turbulence problem itself. It is 
planned to continue the present work along these lines. 

From a more immediately practical viewpoint an under- 
standing of the flow close to a bluff cylinder is important in at 
least two problems, namely, structural vibrations in members 
which themselves shed vortices and structural buffeting 
experienced by members placed in the wakes of bluff bodies. 
Many of these are most appropriately treated by the statis- 
tical methods developed in the theories of turbulence and 
other random processes (ref. 15). These methods are easily 
extended to include the mixed turbulent-periodic phenomena 
associated with problems such as the two mentioned abov.e. 

The research was conducted at GALCIT under the spon- 
sorship and with the financial assistance of the National 
Advisory Committee for Aeronautics, as part of a long-range 
turbulence study directed by Dr. H. W. Liepmann. His 
advice and interest throughout the investigation, as well as 
helpful discussions with Dr. Paco Lagerstrom, are gratefully 
acknowledged. 
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drag coefficient 
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form drag coefficient 
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outside diameter of ring 
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cylinder dimension 

d' 

distance between free vortex layers 

d. 

diameter of ring-supporting wire 
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wake energy 

E\,Ei 

components of wake energy due to periodic 
fluctuations 
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F(n) 

energy spectrum 

Fi{n),F2(n) 

energy spectra of discrete energy 
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continuous energy spectrum 
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tangential velocity in vortex 
Reynolds number 

response characteristic of wave analyzer 
Reynolds number based on ring diameter 
time correlation function 
space correlation function 
distance from vortex center 
radius of vortex 

Strouhal number, based on cylinder dimension, 

UidlUo 

Strouhal number, based on distance between 
free vortex layers, Uid'lUo 
time scale 
time of averaging 
time 

local mean velocity in x-direction 
mean stream velocity 
mean velocity at vortex center 
components of velocity fluctuation 
periodic velocity fluctuations, at frequencies 
7ii and U 2 

random velocity fluctuation 
peak root-mean-square value of velocity fluc- 
tuation 

velocity of vortex relative to fluid 
reference axes and distance from center of 
cylinder 

flatness factor of probability distribution, 

strength (circulation) of a vortex 
Dirac delta function 
positive number 

distance between two points, measured in 
2 -direction 

dimensionless frequency, ^ n. 

U o 

TT 'T 

dimensionless “time” in life of vortex, ^ 

U* d 

dummy variable 
kinematic viscosity 
a value of u 
density 

skewness of probability distribution, 
time interval 

dimensionless spectrum, UoFr(n)IL 
half band width of wave analyzer 

GENERAL CONSIDERATIONS 


Except for the parameters directly related to the sheading 
frequency, the quantities measured were essentially those 
that are standard in turbulence investigations (cf. refs. 12 to 
14). These are briefly reviewed below with some modifica- 
tions required to study the periodic featured. 


REFERENCE AXES 

The origin of axes is taken at the center of the cylinder 
(fig. 1); X is measured downstream in the direction of the 
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Figure 1. — Sketch of GALCIT 20-inch tunnel. 


free-stream velocity, z is measured along the axis of the 
cylinder, which is perpendicular to the free-stream velocity, 
and y is measured in the direction perpendicular to (a:, y ) ; 
that is, y=0 is the center plane of the wake. The free- 
stream velocity is Uo and the local mean velocity in the 
^-direction is U. The fluctuating velocities in the x-, i/-, and 
2 -direction are Uj t;, and tc, respectively. The flow is con- 
sidered to be two dimensional; that is, mean values are the 
same in all planes 2 = Constant. 

SHEDDING FREQUENCY 

The shedding ^ frequency is usually expressed in terms of 
the dimensionless Strouhal number S=UidlUoj where rii is 
the shedding frequency (from one side of the cylinder), d is 
the cylinder diameter, and Uo is the free-stream velocity. 
The Strouhal number. *5 may depend on Reynolds number, 
geometry, free-stream turbulence level, cylinder roughness, 
and so forth. The principal geometrical parameter is the 
cylinder shape (for other than circular cylinders, d is an 
appropriate dimension). However, cylinder-tunnel con- 
figurations must be taken into account, for example, blockage 
and end effects. In water-channel experiments surface 
effects may have an influence. Usually the geometrical 
configuration is fixed, and then S is presented as a function 
of Reynolds number R. 

Instead of Strouhal number it is sometimes convenient to 
use the dimensionless parameter F=nid^lVy where v is the 
kinematic viscosity. 

ENERGY 

The experiments to be described are concerned mainly 
with the velocity fluctuation in the wake, and especially with 
the corresponding energy. 

The energy of the velocity fluctuation at a point in the 

fluid is ^ p(u^ + vi^+y/) per unit volume, where (ti, r, w) is the 

fluctuating velocity and the bar denotes an averaging (sec 
the section “Distril)ution Functions^’). In these experi- 
ments only the component u was measured, and the term 
‘^'nergy^’ is used to denote the energy in tliat component 
only. 

The energy intensity is defined as (u/U„)^. Since the 
mean flow is two dimensional, the intimsity do('s not vary in 
the 2 -direc,tion. At any downstream position in the wake it 
varies in the 7/-direc,tion, normal to the wake. The integral 
of the intensity over a ])lane normal to tlu'. free stn'.am (per 
unit span) is called the wake energy K: 

* Tin* tcTin “shoddinc” is used throughout this n*|>ort, for convt‘nH*n(!o; it is not inoaiit 
to imply anythiiiR about tho mechanism of the formation of fre«^ vortices. 


4 


REPORT 1191 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 






The velocity fluctuation in the wake of a shedding cylinder 
displays a predominant frequency (as well as harmonics) 
which is the shedding frequency. However, except in a 
small Reynolds number range, the fluctuation has random 
irregularities * ^superimposed” on it; that is, it is not purely 
periodic, in the mathematical sense. However, it is con- 
venient to speak of the ‘Aperiodic” and ' ‘random” or turbu- 
lent parts of the fluctuation.® The energy may be written 

U^ = Ut^ + Ui^ + U2^ ( 2 ) 

where Ur^ is that portion of the energy contributed by the 
random (turbulent) fluctuation > is contributed by the 
periodic fluctuation at the shedding frequency rii, and 
corresponds to twice the shedding frequency nt=2ny. (The 
center of the wake feels the influence of vortices from both 
sides and 712 is prominent there, at least near the “beginning” 
of the wake. Higher harmonics are foimd to be negligible.) 

Eqjjation (2) is a kind of spectral resolution, in which 
and U 2 ^ are the energies at the specific frequencies Uy and 712 . 
This type of resolution is called a discrete, or line, spectrum. 
But Ur^ is not a discrete spectral component, for it is the 
energy in the turbulent part of the fluctuation and contains 
“all” frequencies. It has a continuous frequency distribu- 
tion of energy, for which a slightly different definition of 
spectrum is appropriate. This is postponed until the follow- 
ing section. 

Corresponding to equation (2), an equation may be writ- 
ten for the wake energy E and its turbulent and periodic 
components: 

E=Er-\-Ey-\-E2 (3) 

Of particular interest will be the fraction of discrete energy 
{Ey-\-E<^IE at various stages of wake development. 

CORRELATION FUNCTIONS; SPECTRUM 


For r=0, equations (4) and (7) give 

00 

i?t(0)= J F(n)dn=l (8) 

where F{n) is defined as the energy spectrum; that is, F{n) 
dn is the fraction of the energy in the frequency interval n 
to n+dn. It is the fraction of energy “per unit frequency,” 
as contrasted with the discrete energy spectrum discussed in 
the section “Energy.” 

In studies of isotropic turbulence, at Re^molds numbers 
corresponding to those in the present experiments, it is 
found that the energy spectrum is well represented by the 
form 

or, what amounts to the same thing, that the correlation 
fimction is of the form 

7?,(r)=e“^" (10) 

If the normalizing factor K=UolL is used in equation (10), 
L being a characteristic length, then equation (6) gives 

UoF{n) _ 4 . ^ ^ 

L 1+4^2 (LV/C/o") ^ 


which may be conveniently written in terms of the dimen- 
sionless parameters 



v=UoF(,n)lL 

(12a) 

and 

L 




(12b) 

Then 

_ 4 

^ 1 + (2tt/)® 

(11b) 


It is clear from equations (5) and (10) that L is a length 
scale related to the time scale by 


Definitions. — The time correlation function of the fluctu- 
ation u{t) is defined by 






( 4 ) 


where t is a time interval. The time scale is then defined 

by 

7’=J^ R,(It ( 5 ) 


The Fourier transform of Ut defines another fimction 


Then, also 


^ 00 

F(n)=4 J Rt{r) cos 2TUTdT 

r* 00 

Ri{t)= \ F{n) cos 2irnTdT 


( 6 ) 

(7) 


* A turbulent fluctuation is an irrcKular variation, with respect to time, which is charac- 
terized in particular by its randomness and absence of i)oriodicity (cf. ref. 13, p. y). 


L=UoT (13) 

Equation (11b) is used as a. convenient reference curve to 
compare the measurements reported below. 

Periodic functions. — The energy spectrum F(n) is par- 
ticularly well suited to turbulent fluctuations, for which 
the energy is continuously distributed over the fre- 
quencies. For periodic fluctuations the discrete, or line, 
spectrum is more appropriate, but in the present “mixed” 
case it is convenient to wnte the discrete energy, also, in 
terms of F(n). This may be done by using the Dirac delta 
function 5(/i). Thus the energy at the shedding frequency 


7i y IS 




Ui=ui‘^ 5(tt— «i) dn 

(14) 

that is. 




Fi{n) = &{n—ni) 

(15) 


Then the mixed turbulent-periodic fluctuations in the wake of 
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a shedding cylinder are considered to have an energy spec- 
trum which is made up of continuous and discrete parts 
(cf . eq. (2) and appendix A) : 

f F{n) dn=Ur^ f Ft (n) f Fi{n) dn-\- 

Jo Jo Jo 

00 

U 2 ^\ F 2 iji) dn (16) 

that is, 

F{n)=^ FAn)+% 6(n-n,)+% 5{n-n,) (17) 

Space correlation function; phase relations. — The cor- 

relation function defined in equation (4) describes the time 
correlation. Another correlation function which is useful in 
the present study is one which relates the velocity fluctua- 
tions at two points in the wake, situated on a line parallel 
to the cylinder. This is defined by 

( 18 ) 


DISTRIBUTION FUNCTIONS 

Random functions. — The probability density p ($) of 
a random function Ur(t) is defined as the probability of 
finding u, in the interval It may be found by 

taking the average of observations made on a large number 
(ensemble) of samples of Ur{t)y all these observations being 
made at the same time i. This is called an ensemble aver- 
age. If Ur{t) is a stationary process, as in the present case, 
then appeal is made to the ergodic hypothesis and the en- 
semble average is replaced by the time average, obtained by 
making a large number of observations on a single sample 
of Ur(t). The probability density p(() is the number of 
times that Ur is found in (J,J+dJ) divided by the total num- 
ber of observations made. In practice, time averages are 
more convenient than ensemble averages. The averaging 
time Ta must be large enough so that a statistically significant 
number of observations are made. This imposes no hard- 
ship; it is sufficient that Ta be large compared with the time 
scale T, If necessary, the error can be computed. 

Experimentally, ^?U) may be determined by the principle 
illustrated below: 


where f is the distance between the two points. The corre- 
sponding scale is 

(19) 

The function Rz should be particularly suited to studying 
turbulent development. Close to the cylinder it should 
reflect the regularity connected with the periodic shedding, 
especially in a regular, stable vortex street, in which there 
are no turbulent fluctuations. When there are turbulent 
fluctuations and, especially, far downstream where there is 
no more evidence of periodicity, Ez should be typical of a 
turbulent fluid; that is, the correlation should be small for 
large values of f. 

Tho function Rz may be obtained by standard techniques 
applied to the two signals u{Zyt) and u( 2 +f,t). One well- 
known visual method is to apply the signals to the vertical 
and horizontal plates, respectively, of an oscilloscope and to 
observe the resulting ^‘correlation figures^ ^ (or ellipses) on 
the screen (ref. 16). If the signals u{t) are turbulent fluctua- 
tions, then the light spot moves irregularly on the screen, 
forming a light patch which is elliptic in shape. The correla- 
tion function is given by 



71 


^2 


p(i) ~ y 

(21) 

t- 

‘ \dur/dt\ 

(22) 


The most elementary application of this principle is a gi’aphi- 
cal one using a photographic trace of Ur(t). More conven- 
iently, electronic counting apparatus is employed (see the 
section “Statistical Analyzer’'). 

The statistics of Ur{t) are usually described in terms of 
the moments of p(J) and certain functions derived from the 
moments. The moment of order k is defined as 


R. 


a^+h^ 


( 20 ) 




(23) 


where a and h arc the major and minor axes of the ellipse. 

If u{Zyt) is a periodic function, in both time and space, 
then the conelation figure is an elliptical loop (Lissajous 
figure) whose major and minor axes again give Rz according 
to equation (20). Such a case would exist if the wake had 
a spanwise periodic structure. Then Rz(X) would be periodic. 
A special case of this is A'r(f) = l, as would be expected in a 
vortex stn'et, iirovided the vortex filaments are straight and 
parallel to the cylinder and do not “wobble." 


Another useful definition is 

A^»=J^JI*|p(«rf{ (24) 

where Nu is equal to M*. for even values of k. If is 

symmetrical, then A/* is 0 for odd values of k but Nk is not. 

00 

From tlie definition of p{^) it follows that A/o= I p{i;) d^= 
1; J will be normalized by requiring that AA*=l/2, that 








■r : 


is, the mean-square value t/r^=l/2. 

Three useful functions derived from the moments are 



N, N, 

(25) 

Skewness 

M, 

(26) 

Flatness 

Mt 

“ M2^ 

(27) 


The probability density of a function which is partly 
periodic and partly random is expected to display the 
transition from one type to the other. The tendency toward 
the random probability density should be strong. For 
instance, random fluctuations in the amplitude of a sine 
wave result in a large increase in the higher moments. It 
is interesting to study the relation between probability 
functions and spectra, particularly the case where most of the 
energy is discrete but the fluctuation amplitude is random. 

EXPERIMENTAL DATA 


Periodic functions. — The above definitions may be ex- 
tended to the case of a periodic function The 

probability density can be completely determined from a 
single wave length of that is, it is sufficient to take Ta 

equal to the period. This complete a priori information is 
a basic difference between periodic and random functions.^ 

If Ui{t) is measured experimentally then U{^) in equation 
(21) can also be measured. If Ui{t) is given in analytic 
form then U(^) may be calculated from equation (22). Thus 
the distribution densities for simple wave shapes are easily 
calculated. Table I gives the probability densities and 
moments for the triangular wave, sine wave, and square 
wave. Also included is the Gaussian probability density, 
which is a standard reference for random functions. 

The moments of the probability densities of these wave 
shapes are shown in figure 2. The moments for the random 
function increase much faster than those for the periodic 
functions. This results from the fact that the maximum 
values of a periodic function are fixed by its amplitude, while 
for a random function all values are possible. 



Figure 2. — Amplitude distribution moments. 


* For a periodic function the crgodic principle may not be invoked; the ensemble aver- 
aj?e and the time average are not the same (unless the members of the ensemble have random 
phase differences). It is the time average that is computed here, for comparison with the 
experimental results, which are also time averages. 


WIND TUNNEL 

The experiments were all made in the GALCIT 20- by 
20-inch low-turbulence tunnel (fig. 1). The turbulence level 
is about 0.03 percent. The wind velocity’ may be varied 
from about 50 centimeters per second (1 mile per hour) to 
1,200 centimeters per second (25 miles per hour). 

CYLINDERS 

The cylinders used in the experiments varied in diameter 
from 0.0235 to 0.635 centimeter. Music wire or drill rod 
was used. The diameter tolerances are about 0.0002 centi- 
meter. The cylinders spanned the tunnel so that the length 
in all cases was 50 centimeters (20 inches); the cylinders 
passed through the walls and were fastened outside the 
tunnel. 

RINGS 

Some studies were made of the flow behind rings. These 
were made up of wire. Each ring was supported in the 
tunnel by three thin support wires, attached to the ring 
circumference at 120° intervals. Table II gives the dimen- 
sions of the rings used (where d is the wire diameter, D, the 
ring diameter, and d,, the diameter of the support wire). 

VELOCITY MEASUREMENTS 

Velocities higher than about 400 centimeters per second 
were measured with a pitot tube, calibrated against a 
standard. The pressures were read on a precision manom- 
eter to an accuracy of about 0.002 centimeter of alcohol. 
Velocities lower than 400 centimeters per second were 
determined from the shedding frequency of a reference 
cylinder (0.635 centimeter), as explained in the section ^^Use 
of Shedding Frequency for Velocity Measurements.^^ 

Fluctuating velocities were measured with a hot-wire 
anemometer (1/20 mil platinum). Only u{t), the fluctuating 
velocity in the flow direction, has been measured so far. 
The hot-wire was always parallel to the cylinder. 

TRAVERSING MECHANISM 

The hot-wire was mounted on a micrometer head which 
allowed it to be traversed normal to the wake and positioned 
to O.OOl centimeter. The head was mounted on a horizontal 
lead screw which allowed traversing in the flow direction, in 
the center plane of the tunnel. The positioning in this 
direction was accurate to about 0.01 centimeter. The hori- 
zontal lead screw could be turned through 90° to allow 
traversing parallel to the cylinder, for correlation or phase 



THE DEVELOPMENT OP TUEBULENT WAKES FROM VORTEX STREETS 


7 





.■'1 
' *’ 

m 


Vi: 







measurements (see the section ‘^Space correlation function; 
phase relations^O- For this purpose, a second micrometer 
head with hot-wire could be set up in a fLxed position along 
the line of traverse of the first hot-wire. Then correlations 
could be measured between this point and the movable one. 

ELECTRONIC EQUIPMENT 

The hot-wire output was amplified by an amplifier pro- 
vided with compensation up to 10,000 cycles per second. 
The amplifier output could be observed on an oscilloscope 
screen or measured on a Hewlett-Packard Model 400c vac- 
uum-tube voltmeter. Values of were obtained by reading 
the root-mean-square voltage on the voltmeter. (This 
voltmeter is actually an average-reading meter; it reads 
true root-mean-square values only for a sine wave. A few 
of the indicated root-mean-square values, for turbulent 
velocity fluctuations, were checked against true root-mean- 
square values as obtained from the statistical analyzer (see 
the section ^^Statistical Analyzer^^); these may differ up to 
10 percent, depending on the wave shape, but, at present, 
no corrections have been made, since the absolute values 
were not of prime interest.) Usually only relative values of 

were required, but absolute values could be determined 
by comparing the voltage with that obtained by placing the 
hot-wire behind a calibrated grid. 

The frequencies of periodic fluctuations were determined 
by observing Lissajous figures on the oscilloscope; that is, 
the amplifier output was placed on one set of plates and a 
known frequency on the other. This reference frequency 
was taken from a Hewlett-Packard Model 202B audio oscil- 
lator, which supplied a frequency within 2 percent of that 
indicated on the dial. 

FREQUENCY ANALYZER 

Spectra were measured on a Hewlett-Packard Model 300A 
harmonic wave analyzer. This analyzer has an adjustable 
band width from 30 to 145 cycles per second (defined in 
appendix A) and a frequency range from 0 to 16,000 cycles. 
The output was computed directly from readings of the volt- 
meter on the analyzer. It was not felt practicable to read 
output in the frequency range below 40 cycles; therefore, 
the continuous spectrum was extrapolated to zero frequency. 

To determine the discrete spectrum in the presence of a 
continuous background some care was required. In such 
cases the analyzer reading gives the sum of tlie discrete 
spectral energy and a portion of that in the continuous 
spectrum, the proportions being determined by the response 
characteristic of the wave analyzer. The value in the con- 
tinuous part was determined by interpolation between bands 
adjacent to the discrete band and subtracted out to give 
the discrete value, as outlined in more detail in appendix A. 

STATISTICAL ANALYZER 

The statistical analyzer, designed to obtain prol>ability 
functions, operates on the principle described in (he section 


^^Distribution Functions;^^ here u{t) is a voltage signal. A 
pulse train (fig. 3) is modulated by u{t) and is then fed into 
a discriminator which fires” only when the input pulses 
exceed a certain bias setting, that is, only when 
For each such input pulse the discriminator output is a pulse 
of constant amplitude. The pulses from the discriminator 
are counted by a series of electronic decade counters termi- 
nating in a mechanical counter. 

The complete analyzer consists of 10 such discriminator- 
counter channels, each adjusted to count above a different 
value of It will be seen that the probability function 
obtained is the integral of the probability density described 
in the section ‘‘Distribution Functions;” that is, 

P(© = Probability that 14(0 
P(M)d/x 

It is possible to rewrite the moments (see the section “Dis- 
tribution Functions”) in terms of p($), a more convenient 
form for calculation with this analyzer. These are also 
shown in table I. 



Turbulence signal 





Fkjurk 3. — Ropresonlativc signal sociucncc for statistical analyzer. 
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More complete details of the analyzer and computation 
methods may be found in references 17 and 18. 

RESULTS 

SHEDDING FREQUENCY 

Since Strouhal’s first measurements in 1878 (ref. 1) the 
relation between the shedding frequency and the velocity 
has been of interest to many investigators. Rayleigh (ref. 
2, p. 413) pointed out that the parameter UidlUo (now called 
the Strouhal number S) should be a function of the Reynolds 
number. Since then there have been many measurements of 
the relationship (ref. 19, p. 570). One of the latest of these 
is the measurement by Kovasznay (ref. 11), whose determina- 
tion of S(R) covers the range of R from 0 to 10'‘. Kovasznay 
also made detailed investigations of the vortex-street flow 
pattern at low Reynolds numbers. He observed that the 
street is developed only at Reynolds numbers above 40 and 
that it is stable and regular only at Rejmolds numbers 
below about 160. 

The present measurements of S(R) are given in figures 4 and 
5. Except at Reynolds numbers between 150 and 300, the 


scatter is small, and the measurements agree with those of 
Kovasznay. The large number of cylinder sizes used results 
in overlapping ranges of velocity and frequency so that 
errors in their measurement should be ^'smeared’' out. It 
is believed that the best-fit line is accurate to 1 percent. 

The measurements are corrected for tunnel blockage, but 
no attempt is made to account for end effects. With the 
cylinder sizes used no systematic variations were detected. 

NATURE OF VELOCITY FLUCTUATIONS 

It was observed, as in Kovasznay’s work, that a stable, 
regular vortex street is obtained only in the Reynolds number 
range from about 40 to 150. The velocity fluctuations in 
this range, as detected by a hot-wire, are shown on the 
oscillograms in figure 6, for a Reynolds number of 80. These 
were taken at two downstream positions, x/rf=6 and 48, and 
at several values of y/d. (The relative amplitudes are 
correct at each value of x/rf, but the oscillograms for x/d=^S 
are to a larger scale than those for x/rf=6.) The frequencies 
and amplitudes are quite steady; it is quite easy to determine 
the frequencies from Lissajous figures (see the section 
“Electronic Equipment^^, which, of course, are also steady 
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Figure 4. — Strouhal number again.st lUiyriolds fiurnber for circular cylinder. 
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Figure 6. — Oscillograms for R = S0 and ri= 0.1 58 centimeter. 
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Another example, at A^=145, is shown in figure 7(a). 
(The double signal was obtained for correlation studies and 
is referred to later in the section “Spanwise Coi*relation and 
Phase Measurements.” The dotted nature of the trace is 
due to the method of obtaining two signals on one screen, 
using an electronic switch.) 

At Reynolds numbers between about 150 and 500 there 
are irregular hursts in the signal. An example is shown in 
figure 7 (b), at //=180 and xld=i\. 1die hursts and irregu- 
larities become more violent as I\ increases. It is radher 
difficult to determine the frecpiencv. The Lissajous figure 
is unsteady because of the ii*r(‘gularity, hut, in addition, the 
frequency, as wt'll as it can Ix' dctc'rmiiK'd, varies a little. 
This is the reason for the scattc'r in this Reynolds number 
range. 4\v'o sepai*ate plots of S{lk) obtained in two dillerent 
runs are shown in figure 8. They illustrate the erratic 
behavior of S{lk) in this range. 

At Reynolds numh(‘rs above 500, signals like that in figure 
7 (c) were' obtaiiu'd (near the Ix'ginning of the wake). This 
is typical of the velocity (hict nations up to the highest value 
of Ii investigated (about 10,000). There are irregularities, 
but the predominant (slu'dding) freciuency is easy to deter- 
mine from a Lissajous figure. The Lissajo\is figure in this 
case is not a steady loop, as it is at A'=40 to 150, but neither 


is it so capricious as that at /i’= 150 to 500, and the matching 
frequenc}’- is quite easily distinguished from the nearby 
frequencies. 

At j'/(f=48, in this range, all traces of the periodicity have 
disappeared and the fluctuations are typically turbulent. 
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(a) r/= 0.0362 centimerer. 

(b) d = 0.'MS ceutiiiieter. 

Fksure 8. — Plot of S(R) for two single runs. 
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REGULAR AND IRREGULAR VORTEX STREETS 

The above observations show that there are three charac- 
teristic Reynolds number ranges, within the lower end of the 
shedding range. These will be called as follows: 

Stable range 40< 1 50 

Transition range 150</;<300 

Irregular range 300</?< 10,000 + 

As noted above, the actual limits of these ranges are somewhat 
in doubt and may depend on configuration, free-stream tur- 
bulence, and so forth. Also the upper limit of the irregular 
range is undoubtedly higher than 10,000. (Periodic fluctua- 
tions in the wake have been observed up to the critical 
Reynolds number, about 200,000, but the present measure- 
ments did not extend beyond 10,000.) 

In addition to the differences in the nature of the velocity 
fluctuations, the ranges are characterized by the behavior of 
the Strouhal number: In the stable range S{R) is rapidly 
rising, in the irregular range it is essentially constant, and 
in the transition range it is ^ Unstable. 

It will be seen in the further results presented below that 
all phases of the wake development are different in the two 
ranges, stable and irregular, and that they are indeed two 
different regimes of periodic wake phenomena. 

RELATION OF SHEDDING FREQUENCY TO DRAG 

The relation between the Strouhal number S{R) and the 
drag coefficient Cd(R) has often been noted (ref. 19, p. 421). 
Roughly, rising values of S{R) are accompanied by falling 
values of Cd{R) and vice versa. 

The relation to the form drag is even more interesting. The 
total drag of a cylinder is the sum of two contributions: The 
skin friction and tlie normal pressure. At Reynolds numbers 
in the shedding range the skin-friction drag is ‘^dissipated'' 
mainly in the cylinder boundary layer, while the pressure 
drag (or form drag) is dissipated in the wake. It may, 
then, be more significant to relate the shedding frequency 
to the form drag, both of which are separation phenomena. 
The 7+dependence of the pressure drag coefficient 
taken from reference 19, page 425, is shown in figure 5. It 
has several interesting features: 

(a) Cop is practically constant, at the value Cr>p=\. 

(b) The minimum point A is at a value of R close to that 

at which vortex shedding starts. 

(c) Tlie maximum point B is in the transition range. 

(d) In the irregular range Cop{R) is almost a “mirror re- 

llection” of S{R). 

Since the drag (‘oefficient is an “integrated" phenomenoiu 
it is not expected to display so sharply detailed a dependence 
on R as does the Strouhal number, ])ut these analogous varia- 
tions are believed to be closely related to the ])osition of the 
boundary-layer separation point, to which both the shedding 
fre(iuency and the pressure drag are (piite sensitive. 

USE OF shedding; FRE<?UEN<’V for VELOC ITY MEASUREMENTS 

The n'lnarkable dependence* of the shedding frequenev on 
tlu' velocity and the possibility of accurately measuring aS(/?) 


make it possible to determine flow velocities from frequency 
measurements in the wake of a cylinder immersed in the flow. 
At normal velocities the accuracy is as good as that obtain- 
able with a conventional manometer, while at velocities 
below about 400 centimeters per second it is much better. 
(For instance, at a velocity of 50 centimeters per second the 
manometer reading is only about 0.001 centimeter of alco- 
hol.) In fact, in determining S{R) in the present experi- 
ments, this method was used to measure the low velocities by 
measuring the shedding frequency at a second reference c.yl- 
inder of large diameter. The self-consistency of this method 
and the agreement with Kovasznay's results are shown in 
figure 4. 

For velocity measurements it is convenient to plot the 
frequency-velocity relation in terms of the dimensionless 
parameter F (see the section “Shedding Frequency") as has 
been done in figures 9 and 10. The points on these plots 
were taken from the best-fit line in figure 4. They are well 
fitted by straight lines 

(la) 7^=0.212/7-4.5 50</?<150 

(lb) 7^=0.212/7-2.7 300</7<2,000 


which correspond to 

(2a) S=0.212 (1-21.2//7) 50</7<150 

(2b) S=0.212 (1-12.7//7) 300<7?<2,000 


Line (2b) has been plotted in figure 4 to compare with what 
is considered the best-fit line. The agreement is better than 
1 percent. If line (2b) is extended up to /7= 10,000, the maxi- 
mum error, relative to the best-fit line, is 4 percent. 
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Figure 10. — Plot of F againsi R (50<^ 2,000). 

The plot of F(^R) is used as follows: The shedding frequency 
is observed and F=^ri\d?lv is calculated {v is easily determined) ; 
R is found on the F{R) plot and the velocity is calculated 
from R=Uodlv. 

WAKE ENERGY 

From the velocity traces on the oscilloscope (figs. 6 and 7) 
it is clear that in the regular range the fluctuating velocity 
K{t) is purely periodic, while in the irregular range some of 
the fluctuations are random. This difference is illustrated 
in figure 1 1 which shows the distribution of energy intensity 
(ujCof across the wake at two Reynolds numbers, one in the 
regular range, at 77=150, and one in the irregular range, at 
77=500. Onl}^ half the wake is shown for each case; the one 
at 77= 150 is plotted on the left side of the figure and the one 
for 77=500 on the right. 

The total energy intensity {ulUoY at each point was de- 
termined directly from the reading on the root-mean-square 
voltmeter (see the section ^‘Electronic Equipment’^). The 
components at the frequencies rti and 712 were determined 
by passing the signal through the wave analyzer. The curves 
in each half of figure 11 satisfy the equalities 


(0=(t?.)+(k) 




77=150 


77 = 500 


The values of (uI(\Yj (^l/7^)^ and were obtained 

by measurement (and at 77=150 are self-consistent) while 
(Urjl was obtained by difference. The absolute valium in- 
dicated are somwhat in doubt since the vacuum-tvibe volt- 
meter is not a true root-mean-square meter but an* believed 
accurate* to about 10 percent. 

The particular feature illustrated in figure 1 1 (alr(*ady ob- 
vious from the oscillographs) is the absence of turbulent 
en(*rgy at 7i^=150 as contrasted with the early app<*aranc(* of 
turbulent energy at 77=500. This contrast is typical of the* 
r(‘gular and irn*gular ranges. 

The measurements shown were made at 6 diameters down- 
stream, but the same features exist closer to the cvlind(*r. 
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(a) /^=150. 

(b) /? = 500. 

Figure 11. — Wake ejiergy. rf = 0.190 centimeter; xld=Q. 

DOWNSTREAM WAKE DEVELOPMENT 

The downstream development for the case of figure 11 
(but 77=500 only) is sh own in figure 12. The distribution of 
total energA" intensity (u/l^Y is show n on th e left of the 
figure and the discrete energV' intensity (uilUoY, at the shed- 
ding frequency, is shown on the right. Traverses were 
made at 6, 12, 24, and 48 diameters downstream. The 
<liscrete energy decays quite rapidly and is no l onger m easur- 
able at 48 diameters. (Note that the plot of (uJFoY at 24 
diameters is shown magnified 10 times, for clarity.) A plot of 
{U 2 IU 0 Y has not been included since it can no l onger be 
measured at even 12 diameters. The distribution of {Urll'oY 
may be obtained from these curves by difference. 

Figure 13 presents the downstream wake development in 
another way. The w ake energ>^ E was calculated by integra- 
tion of curves like those in figure 12 (cf. the section “En- 
erg> ”) ; that is. 




Figure 13 is a plot of the energy ratio {Ex + E 2 )!E, that is, the 
ratio of the discrete energy relative to the total energy. 

In the irregular range the energies were computed in this 
way at 77=500 and 4,000 (two (*ylinder sizes in each case) 
and 77=2,900 (one cylinder). Figure 13 shows that the 
decay in all these cases is similar and the Avake is completely 
turbident at 40 to 50 diameters. 

The value of x/d for which E^jE becomes zero was deter- 
mined tor a variety of cylinders, varying in size from 0.00 
to 1.3 <*entimeters and at Reynolds numhers from 200 to 
10,000. The value was fo\md to lie between 40 and 50 in all 
cases hut closer to 40. A precise determination is difficult 
(and not important) l)ecause of the asymptotic approach of 
EJE to zero (E 2 is alr(*ady zero at less than 12 diam(*t(*rs). 

In contrast with this, the stal)le range (/f=50 and 100 in 
fig. 13) has no development of turl)ul(*nce Ix'fore 50 diameters. 
1'he |)lots for R= 1 50 and 200 illustrate the rather spectacular 
transition from tin* stabh* rang(* to tin*, im'gular. 


L 
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(a) Total energy intensity {ulUo)^. 

(b) Discrete energy intensity (uiIUq)^. 

Curve for x/d=24 magnified 10 times for clarity. 

Figure 12. — Wake development. d=0,190 centimeter; R = 500. 



For R=50 and 100 the energy ratio remained constant at 
unity up to x/d=l00. Beyond that tlie energy intensity 
is so low that the tunnel turbulence cannot be neglected. 

MEASUREMENTS OF SPECTRUM 

Figure 14 shows spectrum measurements at 0, 12, 24, and 
48 diameters downstream at a Reynolds number of 500, 
in the irregular range. 'Jdie lat(‘ral position y/d chosen for 
the measurement at each x/d is the one for which {ted ( ‘ o)'^ is a 
maximum (cf. fig. 12). The method of plotting is as follows. 
The curve through the experimental points is tin* continuous 
spectrum Ft{u). plotted in normalized coordinat(‘s. 1die 
discrete energies Fi = 6{n—nd and /^2 = 5(n— are indicated 
by narrow ‘d)ands” which should have zero width and infinite 
height but are left '‘open” in the figure. The relative 
energies represented by the areas under the continuous curve 
and under the delta functions, respectively, are marked in 
the figure with values of uditid and wdjtF. 



Figure 14. — Downstream development of spectrum. d= 0.190 
centimeter; /?=500; rii = 440. 


To normalize the continuous spectrum the dimensionless 
parameters ^r{n) and ^ used. In each case 


the curve 


; is included for reference. The nor- 


l + (27r7;)2^ 

malizing coefficient L was determined as follows: 

(a) F’r(O) was found by extrapolation of the measured 
values to n=0. 

(b) F;.(0) and the other values of Fr{n) were normalized 
to make S F{n)dn—\, 

(c) L was found from ^ F,.(0)=4. 


In short, the measured curve and the reference curve were 
made to agree in and in area. This requirement deter- 
mines L. 


In these coordinates the shedding frequency shows an ap- 
parent increase downstream; this is because the normalizing 
parameter L increases. F^or xjd = AS the shedding frequency 
(i. e , Tij) is marked with a dash; it contains no discrete energy 
at this value of xjd. 

The "bumps” in the continuous spectrum, near and 
Uo, indicate a feeding of energy from the discrete to tiu^ 
continuous spectrum. The portion of the spectrum near 
n = 0, which is established early and whi(‘h contains a lai’ge 
part of the turl)ulent energy, seiuns to be unrelated to the 
shedding frecjuency (cf. fig. 15). As the wake develops the 
energy in the bumps is rapidly redistriliu ted (part of it 
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Figure 15. — Downstream development of spectrum, ^/ = 0.953 
centimeter; 7?=4,000: /?i=144. 


decays) to smooth the spectrum, which, in tlie fully de- 
veloped turbulent wake at 48 diameters, tends toward the 

, • • 4 

characteristic curve , "^ 2 * 

1+(27T77)^ 

In figure 16 the spectrum for x!d=\2 and ?//(/=(). 8 is 
plotted together with the one at ///(/= 0. 4^110 curves are sim- 
ilar at low frequencies (large eddies) and at high frecpiencies; 
they differ only in the neighborhood of the discrete band. 
(The slight discrepancy in Ui between this figure and fig. 14 
is due to the fact that they were measured at two different 
times, when the kinematic viscosity v differed. This re- 
sulted in different values of Aq at the same I\.) 

A similar downstream developnu'nt is shown in figure 15 
for /i =4,000. Here the spectrum at jA/ = 6 is smoother than 
that in the previous e.xample (fig. 14). This effect may be 
due not so much to the higher Reynolds number as to the 
fact that the shedding freciuency is closer to the low frequen- 
cies; that is, the shedding freciuency is “embedded” in the 
low-frecpiency turbulent band. It seems to result, at 48 
diameters, in a much closer ai)proach to the ndcTcnce curve. 

Figure 17 shows the spectra at 4S diameters for three 
cvlind(‘rs and several values of //;>/. It is remarkable that 
//-4,000, c/- 0.477 centimeter agrees better with A' = 500, 
({ - 0.100 centimeter than with A’ = 4,000, J = 0.053 ceiiti- 
imOer. This seems to bear out th<‘ above remaj'k about the 
relative influence of R and for the respective shedding 
fn*(iii(*nci(*s arc' 565, 440, and 144. 



Fkjure 15. — Spcrtra at 12 diameters, r/ ^0.190 centimeter: /? = 500; 

a 1 = 4:^0. 


Finally it may be noted that values of which in figure 
1 1 were obtained by difference, cdieck well with the values 
computed from f Vri:n)dn (before normalization of 

Fr{n)). 

Spectra for the regular range are not presented, for they 
are simple discrete spectra. 

SPANWISE CORRELATION AND PHASE MEASUREMENTS 

The function Rz was not measured, but the main features 
of thespanwise correlation* are illustrated in figures 7 and 18. 

h^igure 7 shows thiTe e.xamples, in c'ach of which simulta- 
neous signals were obtainc'd from two hot-wires at xld=i) 
and i//d=\ and separated b\- 50 diameters spanwise. The 
two signals were ohtainc'd simultaneously on the oscilloscope 
screen by means of an electronic switch. This accounts for 
the dottc'd trace's. 

At A’ ==145 (fig. 7 (a)) the c'orrelation is perfect, but there 
is a phase shift. At A'=180 (fig. 7 (b)) the c-omdation is 
still good, but the individual signals occasionally break 

Mil the rcimiiiidor of this .'Section :i di.stiiiction is iiuule between the terms “correlation 
fiinction*’ and “correlation.” The former refers to the function defined in the .section 
“Space correlation function; ph:vse relations,” while the latter is vised in a looser, descrivi- 
tive .«?cnsc. 
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Figure 17. — Spectra at 48 diameters. 



(0) /? = 80; 
100 . 



(b) /? = 220; 

t/d= 10 . 



(c) /? = 220; 

C/d- 105. 



(d) /P = 500; 
C/d 



(e) R-500- 

C/d-\0. 



(f) /?=500; 
C/d - 105. 

L-84891 


Fkiurk 18. — ('orrolation figures. r/ = 0. 158 centimeter; (>; yj(l=\'^ 
exposure, Yz second. 


down. The breakdowns are uiieorrelated at this distance of 
50 dianielers. At /i’=500 (fi^. 7 (e)) (‘Jicli signal still shows 
a predominant frequency. There is some variation in ])hase 
between the two signals. The amplitude irregularities 
appear to be uncorrelated. 


Figure 18 shows the correlation figures obtained by plac- 
ing the signals of the two hot-wires on the horizontal and 
vertical plates, respectively, of the oscilloscope. 

For and f/rf=100 a steady Lissajous figure is 

obtained, showing that the periodic fluctuations at the two 
points (100 diameters apart) are perfectly correlated (but 
they are not in phase). 

For E=220 and 500 there is good correlation only at small 
values of f/rf, that is, only when the two hot-wires are in the 
same ‘'eddy,'' so to speak. For 7^=500 the figures are 
similar to those obtained in fully developed turbulence. 

In obtaining these correlations a remarkable phenomenon 
was observed. The stable vortex street (i. e., 7?<150) has 
a periodic spanwise structure. This was shown by a phase 
shift on the Lissajous figure, as the movable hot-wire was 
traversed parallel to the cylinder. From the phase coinci- 
dences observed, the wave length parallel to the cylinder 
was about 18 diameters at a Reynolds number of 80. It has 
not been determined whether this periodicity structure is 
due to a “waviness" in the vortex filaments or whether the 
vortex filaments are straight but inclined to the cylinder 
axis. 

STATISTICAL MEASUREMENTS 

A few amplitude distribution functions were measured 
and are shown in figure 19. One measurement is in the 
stable range; the other shows downstream development in 
the irregular range. 



KuaTRE 19. — Distribution functions. ^=0.190 centimeter. 


The table in figure 19 shows values of c and a computed 
from these curves. The behavior, of course, is as expected, 
but tlie numerical values are of some interest. These values 
(and the curves) show that at A'=100 the signal was prac- 
tically triangular but had rounded “tops." At 71^=500 the 
downstream development of randomness is shown by tlic 
tendency of c and a toward the Gaussian values. 

The distribution is in fact not Gaussian, as may be seen 
in the figure, for its skewness a is quite high. 

VORTEX RINGS 

The flow behind wire rings was br.!efly investigated. The 
dimensions of the rings used are given in table II. 
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With the rings of diameter ratio Djd=\^ vortices are shed 
from the wire in almost the same way as from the straight 
wire, and there is apparently an annular vortex street for 
some distance downstream. The Stroiilial number, meas- 
ured from 7?=70 to 500, is lower than that for the straiglit 
wire (about 3 percent at 7? ==500 and 6 percent at 7i’=100). 

Fluctuating velocity amplitudes were measured in the 
wake at several downstream positions. The results for the 
largest ring, measured along a diameter, are shown in figure 
20. It should be noted that -yju^ rather than the energy 
has been plotted here (cf. fig. 11); only relative values were 
computed. Close behind the cylinder the wake behind 
the wire on each side of the ring is similar to that behind the 
straight wire, but the inside peaks are lower than the outside 
peaks. This may be partly due to the interference of the 
hot-wire probe, for a similar effect, much less pronounced, 
was noticed in the measurements behind a straight wire. 

Farther downstream there was some indication of strong 
interaction between the vortices, for a peak could not be 
followed ‘‘smoothly’’ downstream. However, the investi- 
gations were not continued far enough to reach conclusive 
results. At about 40 diameters downstream the flow became 
unstable. 



Figure 20. — Shedding from a ring. </=0.1()8 centimeter; D=1.59 
centimeters; 7^=100; wi = 84. 


The ring with 77^/= 5 behaved somewhat difrerently. 
The observed frequencies gave values of Strouhal number as 
shown in table 111. TIk' table shows values of S and 11 
bas(‘d on the wire diameter, as well as valu(*s of Sp and I\r> 
based on ring diamc'ter. B(*tw(*(‘ii 74^=153 and 182 tlu're is 
a su(l(l(‘u inen^ase in S, and at higher Rtwiiolds numb(‘rs, in 
what corresponds to the irregular range, the shedding is 
similar to that from a straight wire, while in the stable range 
the shedding is at a much lowtu* freqiKuicy. From the 
observations made it s(‘ems lik(‘ly that in tlu^ stable' range the 
ring acts like a disk, shedding the vortex loops observed by 
Stanton and Marshall (r(*f. 18, ]). 578, and r(‘f. 20). Stanton 
and Marshall do not give their frequency-vc'locjty observa- 
tions except at the critical Ro, where sluidding first starts. 


They observed this to be at about 7?/)= 200, with a corre- 
sponding Sd of 0.12. 

Again, these experiments were too incomplete to warrant 
definite conclusions, but the difference in behavior for Dld= 10 
and Dld—5 is interesting. This behavior is similar to that 
observed by Spivack (ref. 21) in his investigation of the 
frequencies in the wake of a pair of cylinders which were 
separated, normal to the flow, by a gap. He found that 
when the gap was just smaller than 1 diameter instability 
occurred. For larger gaps the cylinders behaved like indi- 
vidual bodies, while for smaller gaps the main frequencies 
were, roughly, those corresponding to a single bluff body of 
dimension equal to that of the combined pair, including the 
gap. 

DISCUSSION 

The most significant results of this investigation may be 
discussed in terms of the Reynolds number ranges defined 
in the section “Regular and Irregular Vortex Streets,” 
namely, the stable range from 7?=40 to 150, the transition 
range from 77=150 to 300, and the irregular range above 
77=300. 

STABILITY 

The transition range from 77=150 to 300 displays the 
characteristics of a laminar-turbulent transition, and it is 
instructive to compare the stability of the flow around the 
cylinder with boundary-layer stability. The flow in the 
irregular range has turbulent characteristics, while in the 
stable range it is essentially viscous. 

The Reynolds number regimes may be described as follows: 
Below 77 = 40 the flow around the cylinder is a symmetric, 
viscous configuration, with a pair of standing vortices be- 
hind the cylinder. At about 77=40 this symmetric configu- 
ration becomes unstable. It changes to a new, stable con- 
figuration which consists of alternate periodic breaking away 
of the vortices and formation of a regular vortex street. The 
change at 77=40 is not a laminar-turbulent instability; 
it divides tw^o different ranges of stable, viscous flow. In 
either range, disturbances to the stable configuration will 
be (lamped out. 

On the other hand, the transition range from 77= 150 to 300 
involves a laminar-turbulent transition. To understand 
how this transition is related to the vortex shedding, it is 
iK'cessary to know something about the formation of the 
vortic.('s. Involved in this formation is the circulating 
motion behind the cylinder as shown in the following sketch. 
A free vortex layer (the separated boundary layer) springs 
from ('acli se])aratiou point on the cylinder. This five layer 
and the backflow b(‘hind the cylinder ('stahlish a circulation 
from which fluid “breaks away” at ivgiilar intervals. 


Sheer 



Skelcli 2. 
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The laminar-turbulent transition is believed to occur 
always in the free vortex layer; that is, the circulating fluid 
becomes turbulent before it breaks away. Then each vortex 
passing downstream is composed of turbulent fluid. 

The point in the free vortex layer at which the transition 
occurs will depend on the Reynolds number. This transition 
was actually observed by Schiller and Linke (ref. 18, p. 555, 
and ref. 22) whose measurements were made at cylinder 
Reynolds numbers from 3,500 to 8,500. The distance to the 
transition point, measured from the separation point, de- 
creased from 1.4 diameters to 0.7 diameter, and for a given 
Re^molds number these distances decreased when the free- 
stream turbulence was increased. Dryden (ref. 23) observed 
that at some value of /?, depending on free-stream turbulence 
and so forth, the transition point in the layer actually reaches 
the separation point on the cylinder. Transition then re- 
mains fixed and vortex shedding continues, essentially 
unchanged, up to Reynolds numbers above 100,000, that is, 
up to the value of R for which transition begins in the cylinder 
boundary layer ahead of the separation point. It is quite 
likely that even above this critical value of R the phenome- 
non is essentiall}’’ unchanged, but now the vortex layers are 
much nearer together and the vortices are diffused in a much 
shorter downstream distance. 

In summary, vortex formation in the stable range occurs 
without laminar-turbulent transition. The circulating fluid 
breaks away periodicallj^ and alternately from the two sides, 
forming free ^Viscous’’ vortices which move downstream 
and arrange themselves in the familiar vortex street. In 
the irregular range transition occurs in the circulating fluid 
before it breaks away, and the vortices are composed of 
turbulent fluid. The transition range corresponds to the 
similar range in boundary-layer stability, and it displays a 
similar intermittency. The values 7?=150 and 300 used to 
define the range are expected to be different in other experi- 
ments, depending on wind-tunnel turbulence, cylinder 
roughness, and so forth. 

SHEDDING FREQUENCY 

The Strouhal number and Reynolds number dependence 
is different in the two ranges. In the stable range S{R) is 
rapidly rising, while in the irregular range it is practically 
constant. 

Fage and Johansen, who investigated the structure of the 
free vortex layers springing from the separation points on 
various bluff cylinders (ref. 9), made an interesting observa- 
tion on the relation of the shedding frequency to the distance 
between the vortex layers. This distance increases as the 
cylinder becomes more bluff, while the slu'dding frequency 
decreases. In fact, if a new Strouhal number aS" is defined in 
terms of the distance (V between the fn‘e vortex lay(*rs 
(instead of the cylinder dimension d), tln'ii a univiu-sal 
value S"«0.28 is obtained for a variety of (bluff) cyliruh'r 
shapes. The measurements of reference 9 were madi' at 
7^=20,000, but it is believed that the similarity exists ov(*r 
the whole irregular range. It does not extend to the stable 
range. To check this point the shedding frequency was 
measured in the wake of a half cylinder placed with the flat 
face broadside to the flow. It was found that S{R) was 


rising for Reynolds numbers below 300 and then became 
practically constant at the value aS=0.140. For a similar 
case, at 7/=20,000, Fage and Johansen found aS=0.143. 

The universality of the constant S' is useful in systematiz- 
ing the shedding phenomena (at least in the irregular range). 
It indicates that when the circulating fluid behind the cyl- 
inder is turbulent then the formation of free vortices is 
similar for a variety of bluff shapes and over a wide range 
of Reynolds numbers. 

Finally, the relation between Strouhal number and form 
drag coefficient has been mentioned in the section ‘‘Relation 
of Shedding Frequency to Drag.^^ In the irregular range 
the slight variations in S{R) reflect slight variations of Cop 
and so, probably, of the separation point. However, con- 
stancy of Cdp is not enough to insure a fixed separation point. 
For instance, Cop remains practically constant down to 
Reynolds numbers below the shedding range, but the separa- 
tion point there is farther back than it is at higher Reynolds 
numbers. It would seem worth while, and fairly easy, to 
measure the position of separation as a function of Reynolds 
number over the whole shedding range, that is, to complete 
the data available in the literature. 

DOWNSTREAM DEVELOPMENT 

The way in which the wake develops downstream is quite 
different in the stable and irregular ranges. 

When the circulating fluid breaks away before the occur- 
rence of transition in the free vortex layers (i. e., below 
7^=150), then the free vortices which are formed are the 
typical viscous vortices. JJiere is no further possibility 
for the fluid in them to become turbulent. The vortices 
simply decay by viscous diffusion as they move downstream 
(see the section “Spread of Vortex Street’’ in appendix B). 

When turbulent transition does occur, then the vortices 
which are formed consist of turbulent fluid. They diffuse 
rapidly as they move downstream and are soon obliterated, 
so that no evidence of the shedding frequency remains. This 
development to a completely turbulent wake takes place in 
less than 50 diameters. In terms of the decay of the discrete 
energy (fig. 13), the development is roughly the same for 
Reynolds numbers from 300 to 10,000. This again indicates 
a remai'kable similarity over the whole irregular range. 

The stable and irregular ranges are also cliaracderized by 
the diflerence in the (uiergy s])ectra of the velocity fluctua- 
tions. It has be(*n ])oint(Hl out that in the irregular range a 
continuous, or turbuh'iit, ])art of the s])ectrum is establislied 
at the beginning of tlu‘ wake develo])ment.. This turbulence 
is a r(‘sult of th(‘ transition in thc^ fi('(» vortex layers and might 
be ex])ect(‘d to b(‘ in(le])(‘nd('nt (at first) of the periodic part 
of the fillet nation, which results from t he ])eriodie, shedding. 
Ind(‘('d, most of the (UK'/gy at first is (T)neentrated at the 
shedding frecjiKMiey rii (some at // 2 ), and it may be re])resented 
as a diser('t(* (dc'lta function) ]>art^ of the s])(M*trum, within the 
aecnraey of tin* m(‘asm(*m(*nls (ef. ap])(*ndix A). IIowi'vit, 
th(* continuous and diser(‘tc‘ ])arts an' not ('iitirely inde])end- 
('iit, as shown by the bum])s lU'.ar ui and ri 2 (fig. 14). This 
may he regarded as a result of energy “fei'ding” from the 
disen'ti' to the eonlinuoiis ]>arts of the six'etriim, and it 
])roe('(‘ds in a way which tc'iids to smooth the spi'edrum. Such 
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transfer of energy between spectral bands is a process de- 
pending on the nonlinear terms of the equations of motion. 
The ^^activity^^ in the spectrum, at any stage of its develop- 
ment, may be regarded as an equilibrium between the 
nonlinear and the viscous terms. It is an important problem 
in the theory of isotropic turbulence. 

The spectral activity near the frequency of discrete energy 
might be looked upon as a simplified case in which a single 
band has an excess of energy and the spectral energy flow is 
unidirectional, that is, out of it into the adjacent bands. 
However, the nonhomogeneous character of the field involved 
(the wake) reduces the simplicity, for it is necessary to take 
account of energy transfer across the wake. One interesting 
possibility is to superimpose a homogeneous (isotropic) 
turbulent field, by means of a screen ahead of the shedding 
cylinder, and to study the effect of this field on the spectral 
activity near the discrete band. Although the wake will still 
introduce nonhomogeneity (not even counting the periodic 
part of the motion), it may be possible to arrange the relative 
magnitudes to give significant results from the simplified 
model. 

To study such problems the technique for measuring the 
spectrum (appendix A) near the frequency of discrete energy 
will be improved. 

To summarize, it is suggested that the initial development 
of the spectrum might be regarded as follows: The continu- 
ous and the discrete parts are established independently, the 
one by the transition in the vortex layers and the other by 
the periodic shedding. The turbulence due to the transition 
is the ''primary'' turbulent field and its spectrum is the 
typical, continuous (turbulent) spectrum. (It has been noted 
in the section "Measurements of Spectrum" that the low- 
frequency end of the spectrum is established early; it would 
contain only energy of the primary field.®) The discrete part 
of the spectrum is embedded in the turbulent part, and it 
thereby is "excited" into spectral transfer. Some of its 
energy is transferred to the adjacent frequency bands result- 
ing, initially, in the development of bumps in the continuous 
spectrum. Subsequently, as the spectral transfer proceeds, 
the spectrum becomes smooth. 

The above discussion is an abstract way of saying that the 
vortices are diffused by a turbulent fluid (instead of a viscous 
one). The diffusion involves the nonlinear processes typical 
of turbulence; the study of these processes, in terms of 
spectrum, is an important problem. 

There is a similar case of turbulent, periodic structure in 
(he flow field between two cylinders, one of which rotates. 
Taylor’s discovery of the periodic structure of the flow is 
well known (ref. 24). When the inner cylinder rotates, it 
is possible to obtain a steady, regular arrangement of ring 
vortices, enclosing the inner cylinder, and having, alternately, 
opposite directions of circulation. Above a critical value of 


the speed of rotation this laminar, periodic structure btr- 
comes unstable and the fluid becomes turbulent, but alternate 
ring-shaped vortices still exist at speeds several hundred 
times the critical speed (ref. 25). 

STATISTICS 

The probability distribution functions (fig. 19) display the 
characteristics which are expected, from the other observa- 
tions. The contrast between the functions at /?=100 and 
500, that is, in the stable and irregular ranges, respectively, 
is quite evident. In the irregular range, even at ild=Q, 
where most of the energy is discrete, there is a marked 
irregularity in the fluctuation, as shown by the high value of a. 

However, these descriptions are little better than qualita- 
tive, and it is hoped to obtain more interesting results by 
extending these statistical methods. Of particular interest 
in the development of random from periodic motion would 
be the relation between the probability distributions and the 
the spectra. For instance, it is plain that a purely periodic 
function (discrete spectrum) will have a probability dis- 
tribution with finite cutoff, while development of random 
irregularities in the function’s amplitude is strongly re- 
flected in (1) a "spreading" of the distribution function to 
higher values of f and (2) the appearance of a continuous 
spectrum. However, the relation between the two is not 
unique; that is, the spectrum does not give (complete) 
information about the probability distribution, and vice 
versa. 

SUGGESTIONS FOR FUTURE INVESTIGATIONS 

Some further lines of investigation indicated by these 
experiments are summarized below. 

(a) The transition from the stable to the irregular range 
should be investigated with controlled disturbances, for 
example, cylinder roughness and free-stream turbulence. 
It is expected that the limits of the transition range (roughly 
/^=150 to 300 for the experimental conditions here) will be 
lower for higher free-stream turbulence or cylinder roughness. 
The critical cylinder Reynolds numbers should be related to 
corresponding numbers for the transition point in the free 
vortex layers (based on distance from separation point or on 
the thickness of the layer). 

Such studies of stability to different disturbance amplitudes 
and frequencies are well known in the case of the boundary 
layer. A variation of the experiments of Schubauer and 
Skramstad (ref. 2G), who used an oscillating wire in the 
boundary layer to produce disturbances of definite fre- 
(juencies, would be to use a second shedding cylinder. 

(b) A st\uly of the spectral development in the neighbor- 
hood of a discrete band, the effect of a turbulent field on its 
activity, and so forth (discussed in the section "Downstream 
Development") may be the most* fruitful continuation of 
these experiments. So far, the problem has been approached 
only in the theory of isotropic turbulence, where it has not 
advanced much beyond tlie similarity considerations of 


« In the theM)ry of honioKcnooiis tiirhiiloiico it i.s shown Hint the low-froqurncy riul of 
tho siM'otniin is invjirinnt, a proi-M'rty related to the liOitsianski invariaiil . 
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Kolmogoroff, and very little is known about the form of the 
spectral transfer function. 

Interactions between discrete bands, for example, at 
slightly different frequencies, can be studied by the use of 
two or more cylinders arranged to ‘^interfere’’ with each 
other (some such studies have been made by Spivack (ref. 
21) but not from this viewpoint), or possibly by using one 
cylinder having diameter changes along its span. 

(c) Townsend has recentl}^ used the concepts of intermit- 
tently turbulent flow and local isotropy in his investigations 
of the turbulent wake and has obtained a new description of 
its structure (ref. 27) . His studies were made at downstream 
distances of 80 diameters or more, so that the wake was fully 
turbulent. Probably the structure he describes is essentially 
the same up to the beginning of the fully developed wake 
(about 50 diameters), but then there is the question of how^ 
it is related to the earlier developments. The most obvious 
‘‘early developments’^ are the turbulent transition in the 
free vortex layers and the periodic shedding. (Although the 
shedding frequency is no longer distinguished far down- 
stream, it is prominent in the early spectral developments 
and thus has an influence on the downstream wake.) 

Such studies will involve considerably more detailed 
investigations of the wake structure than were made here, 
possibly along the lines of Townsend’s experiments and the 
classical measurements of energy balance across the wake. 
The other tw^o components of the energy and will be 
needed. 

(d) The nature of the circulating flow behind the cylinder 
and the formation of free vortices, that is, the shedding 
mechanism, should receive further attention. 

(e) The spanwise periodic structure of the vortex street 
should be investigated, beyond the very cursory observations 
made here. In particular, a study of the stability of single 
vortex filaments seems important. 

(f) Measurements of the fluctuating forces on the cylinder, 
due to the shedding, would be interesting and should have 
immediate practical applications. There seems to be very 
little information about the magnitude of these forces. It 


might be obtained either by direct measurement of forces 
(on a segment) or pressures (with pressure pickups) or 
inferred from measurements of the velocity fluctuations close 
to the cylinder. In addition to the magnitude of the force 
or pressure fluctuations, their spanwise correlation is of 
prime importance. 

CONCLUSIONS 

An experimental investigation of the \vake developed 
behind circular cylinders at Reynolds numbers from 40 to 
10,000 indicated the following conclusions: 

1. Periodic wake phenomena behind bluff cylinders may 
be classified into two distinct Refolds number ranges 
(joined by a transition range). For a circular cylinder these 
are: 

Stable range 40</i’<150 

Transition range 150</?<300 

Irregular range 300<//< 10,000 + 

In the stable range the classical, stable Karman streets are 
formed; in the irregular range the periodic shedding is 
accompanied by irregular, or turbulent, velocity fluctuations. 

2. The irregular velocity fluctuation is initiated by a 
laminar-turbulent transition in the free vortex layers w4iich 
spring from the separation points on the cylinder. The 
first turbulent bursts occur in the transition range defined 
above. 

3. In the stable range the free vortices, which move 
downstream, decay by viscous diffusion, and no turbulent 
motion is developed. In the irregular range the free vortices 
contain turbulent fluid and diffuse faster; the wake becomes 
fully turbulent in 40 to 50 diameters. 

4. A velocity meter based on the relation between velocity 
and shedding frequency is practical. 

5. In the stable range a spanwise periodic structure of the 
vortex street has been observed. 

6. An annular vortex-street structure has been observed 
behind rings having a diameter ratio as low as 10- 

California Institute of Technology, 

Pasadena, Calif., May 29, 1052. 


APPENDIX A 


experimental analysis of SPECTRUM 


These notes supplement the brief descriptions in the 
sections “Frequency Analyzer” and “Measurements of 
Spectrum.” 

ANALYZER RESPONSE 

Consider the response of a spectrum analyzer, such as 
that used in the ])resent experiments, to a mixed periodic- 
random input, and in particidar considcT the ])roblem of 
inferring the input from the output. 

'Hie in])u t, an energy or ])ower, has a random and a ])criodic 
component: 


U^ = Ur^-\-Ui^ 


The corresponding spectra arc defined by 



(Al) 


(A2) 
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where 

F,{n) = 5{n-n,) (A3) 

and 8{n) is the Dirac delta function. 

The response characteristic of the analyzer may be ob- 
tained by considering the effect of a periodic input. When 
the analyzer setting coincides with the input frequency Ui 
the output is a maximum, and when the setting is moved 
away from Ui the output falls off. The response character- 
istic is 


acteristic is quite sharp, relative to the frequency intervals 
of interest. The values 30 to 145 seem quite high, but they 
are a little misleading because of the high attenuation used 
to define oj. For example, if the response-characteristic half 
band width w is 30 cycles per second, it has a total width of 
only 6 cycles per second at 50-percent attenuation. 

SEPARATION OF DISCRETE ENERGY 

To separate the discrete energy Ui^ from the continuous 
spectrum the following procedure is used. 


at setting (^4) 

^ ^ ' Output at setting n^=rii 

The output spectrum G{ua) of the analyzer is related to the 
input spectrum F{n) by (cf. ref. 15) 

G(nA)=\ F(n)R(n—nA)dn 

r Fr(n)R(n-nA)dn+'^ ( 6(ri,,-rj^)7?(n-nx)rfn 
tt* Jo Jo 


r Fr{n)R{n-nA)dn+%R{n, 

Jo 


— TIa) 


(A5) 


Since R(n—riA) is sharp, that is, almost a delta function (see 
the section “Half Band Width”), Fr{n) may be considered 
to be constant over the significant interval of integration in 
equation (A5) . Then 


G(nA)=^= FMa) R («,-n J 




u 




where 


Q—^ R{n—'nA)dn==^ R{n— 


riA) d%A 


(A6) 

(A7) 


is the area under the response characteristic. 

Equation (A6) gives the output for a mixed periodic- 
random input. It is required to find the separate terms 
which make up this sum. The procedure is outlined in the 
section “Separation of Discrete Energy” below. 

HALF BAND WIDTH 

The resolution of the analyzer is determined by its half 
band width co. This is defined as the number of “cycles off 
resonance” at which the output falls off to 0.01 percent; that 
is 


/f(n,-w) =0.0001 


(A8) 


For an ideal analyzer the response characteristic would be a 
delta function, but oven with half band widths from .30 to 
145 (which is the jange of the analyzer used here) the char- 


u^G{r)A)0 





Sketch 3. 

At 7 ii + co and rii— co (sec sketch) the contribution from W 
is only 0.01 percent, so the measured points there are assumed 
to lie on the continuous spectrum. It is assumed at first that 
the continuous spectrum between those points may be deter- 
mined by interpolation, and its value at Ui is calculated 
Then u? is determined by difference and the last term in 
equation (A7) is calculated, since the form R{n) is known. 
The first term in equation (A6) then gives the values of 
G{n^ in the vicinity of n \ ; these should check the measured 
values. 

If, however, the continuous spectrum within the band 
width has a bump, then the above calculation is not self- 
consistent, and the true values can be determined by 
successive estimates of 

In principle the method is satisfactory, but in practice the 
accuracy is low because in the regions of interest, that is, 
near peak frequencies, it depends on the differences of 
relatively large quantities. One of these, R(n), is known 
precisely, but the precision is difficult to realize since the 
settings on the analyzer cannot be read accurately enough. 
For the spectral investigations discussed in the section 
‘'Downstream Development” the technique will be improved, 
by monitoring the analyzer with a counter. 


APPENDIX B 


NOTES ON VORTEX-STREET GEOMETRY AND SHEDDING FREQUENCY 


The regularity of the vortex shedding and its sensitivity 
to velocity changes have undoubtedly intrigued everyone 
who has investigated the flow past blpff bodies. However, 
as K^rmdn pointed out in his first papers on the vortex street, 


the problem is inherently difficult, involving as it does the 
separation of the boundary layer from the cylinder, and 
there is yet no adeciuate theoretical treatment of the mech- 
anism. 
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Tlie following notes may be useful as a summary of the 
interesting features of the problem. They are based largcly 
on the literature but include some results obtained during 
the present experiments. Chapter XIII of reference 19 has 
a very useful review and list of references. 


IDEALIZED KARmAN STREET 


Karnnin’s theory treats a double row of potential 
vortices, infinite in both directions. The distance between 
the rows h and the spacing of the vortices in each row I arc 


constants. The vortices have strength (circulation) T which, 
with the geometry, determines the velocity T' of the street 
relative to the fluid. The theory shows that the configura- 
tion is stable when the rows are staggered by a half wave 
length and the spacing ratio is 


=0.281 


(Bl) 


The circulation and velocity relative to the fluid arc then 
related by 



Two of the parameters (A, I, r, and V) must be determined 
from some other considerations. In the real vortex street 
they must be related to the conditions at the cylinder. 

REAL VORTEX STREET 

The real vortex street, even in the stable range, diffei-s 
from the idealized one in the following points: 

(1) The street is not infinite. It starts shortly down- 
stream of the cylimler and eventually loses its identity far 
downstream. However, the classical vortex-street patterns 
extending for 10 or more wave lengths should be a good 
approximation. 

(2) The vortex spacing is not constant. In particular, 
the lateral spacing h increases downstream. 

(:i) The real vortices must have cores of finite radius. 
These grow downstream, so that the vortices diffuse into 
each other and decrease their circulation. For the, same 
reason the velocity V is expected to differ considerably 
from the theoretical value, since it is strongly dependent 
on the configuration. 

Related to these considerations is the way in which the. 
vortices are first formed. At Reynolds numbers below the 
shedding range a symmetrical pair of eddies is formed at 
the ba<-k of the cylinder. As the Reynolds number increases 
these two eddies grow and become more ami more elongated 
in the flow direction, until the configuration is no longer 
stable and become asymmetric. Once, this occurs the 
circulating fluid breaks away ’’ alternately from eacli side 
to form free vortices which flow downstream and arrange 
themselves into the regular, stable vortex street. 


In the irregular range the process is similar, except that 
the fluid is turbulent (because of the transition in the free 
vortex layers). 


DOWNSTREAM VORTEX SPACING 

In the flow past a stationary cylinder the frequency with 
which vortices of one row pass any point is given by 


ai 


(B3) 


This must be the same as the shedding frequency 


Ui=SUold 

Two useful expressions result: 


or 



(B4) 

(B5) 

(B6) 


In a real vortex street, 


T"-^0 far downstream and then 4 

a 



Or, if Sljd is known from measurements, then VlUo may be 
eomputed. 

An example of measured values of Ijd is shown in figure 21. 
These were taken from the streamline plot obtained by 
Kovasznay (ref. 11) at /f=53 (for which *5=0.128). There 
is a little scatter, but Z/rf does approach the constant value 
l/*5=7.8. 


The scatter, while relatively unimportant in the case of 
//d, gives very low accuracy for values of VlUo calculated 
from equation (B6). These have also been plotted in figure 
21. It is surprising that some of the values, near the 
cylinder, are negative (corresponding to values of Ijd higher 
than l/*5); it is believed that this results from the combined 
difficulty of estimating the vortex centers, especially near the 
cylinder, and the sensitivity of equation (B6). (However, 
it must he noted that negative values of V are not impossible. 
Negative V simply means that the vortex velocity is directed 
upstream relative to the fluid, while it is still downstream 
relative to the cylinder. Such a possibility exists at low 
values of x/d, where the mean velocity at the edges of the 
wake is considerably higher than Uo.) 

Another way to obtain VlUo is to assume that the vortex 
centers move with the local mean velocity. Kovasznay's 
paper includes measurements of mean velocity profiles. 
From his results the mean velocity along the line of vortex 


centers 


f/* has been determined and 


from 


It 


V_ 

IK 




U* 

Uo 


has been calculated. The result is plotted in figure 21. 
Near the cylinder it does not agree with the values obtained 
by the previous method; it is believed that this is principally 
due to the difficulties mentioned above and that the deter- 
mination of VlUo from \ — (U*l Uo) is more accurate. 


7 I^ossibly the broakine: away should be reRarded as primary. rosultinR in asymmetry. 
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Uo 


x/d 



Figure 21. — Vortex-street geometry. Calculated from data in 
reference 11. 


LATERAL SPACING 

The lateral spacing, at least initially, must be determined 
by conditions near the cylinder. Tlie way in which tl’is 
spacing increases downstream is discussed, for the stable 
range, in the section ^^Spread of Vortex Street.^’ 

In the irregular range, the dependence of the shedding 
frequency on the distance between the free vortex layers, 
noted by Fage and Johansen (see the section ^^Shedding 
Frequency”), leads to an interesting estimate of the initial 
lateral spacing of the free vortices. The maximum distance 
d' between the free vortex layers, instead of the cylinder 
dimension d, may be used to define a new Strouhal number 

S'=nS (B7) 

C. ' o 

F^age and Johansen found that , whei’cas S varies considerably 
with cylinder shape, iS' is jiractically constant for a variety 
of bluff cylinders. Now the initial lateral spacing Id i>f tl’.e 
free voi'tici's will be roughly the sanu' as r/', possibly a little 
smaller: 


Then, comparing with equations (B6) and (B7) 

1 — € 

From the measurements of Fage and Johansen, S'«0.28. 

The factor - b Thus equation (B9) gives h'll = 

0.28; that is, the spacing ratio agrees with Karman's value, 
at least close to the cylinder. 

SHEDDING FREQUENCY 

There is yet no adequate theory of the periodic vortex 
shedding, and it is not clear what is the principal mechanism 
which determines the frequency. 

The downstream spacing ratio is related to the shedding 
frequency by equation (B3) and to the lateral spacing by a 
stability criterion (e. g., Karman’s value of 0.28 for the 
idealized street). It might be considered that the shedding 
frequency is determined b}" the spacing requirement, or, 
conversely, that the shedding is primary and determines 
the downstream spacing. The latter viewpoint seems the 
more plausible one; that is, the shedding frequency is estab- 
lished by a mechanism which depends on features other than 
the vortex spacing. It is necessary to obtain a better under- 
standing of the flow field near the cylinder. One of the 
elements involves the problem of separation, particularly 
the nonstationary problem. Another that requires more 
study is the flow field directly behind the cylinder. 

With a better knowledge of these, and possibly other, 
features it may be possible to set up a model of the shedding 
mechanism. In the meantime it is not clear whether the 
vortex spacing requirement is decisive in determining the 
frequency. 

DESTABILIZATIONS OF SHEDDING 

The following experiment illustrates the dependence of 
the periodic shedding on “communication” between the free 
vortex layers, that is, on the flow field directly behind the 
cylinder. A thin flat plate was mounted behind the cyl- 
inder in the center plane of the wake (fig. 22). It was 
completely effective in stopping the periodic shedding. 
Spectrum measurements in the flow on one side of the jilate 
are shown in figure 22. At 7?=7,500 no significant fre- 
quencies could be separated out from the continuous back- 
ground. At 7? = 3,200 there were several predominant fre- 
quencies (all higher than the shedding frequency for the 
cylinder), but, by the time the flow reached the end of the 
plate, 5 diameters downstream, it was comi)letely turbulent. 
(The shedding frequency for the cylinder is marked in 
the figures.) 

The important effect, on the shedding, of the flow field 
directly behind the cylinder is a])parent. Probably an even 
shorter length of plate would be etfective in dcstablizing 
the periodic shedding, and there may be a most etfective 


(B9) 


(B8) 


' Thr stabilily consKUMr*! in this sod ion is not with rosin'd to V.unimv-tnvbnlont transi- 
tion; it conoorns tho st:il>ilify of tlu* poriodio shodding (of. tho soction “Stability"). 
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(a) 72 = 3,200; 2. 

(b) 72 = 7,500; x/rf=l. 

Figure 22. — Effect of downstream plate on wake frequencies. 


position for such an interference element. Kovasznay 
remarks that the hot-wire probe used in investigating the 
vortex street must be inserted from the side, for if it lies in 
the plane of the street it has a strong destablizing influence. 

A more complete study of the destabilization of shedding 
by such interference devices may be quite useful from a 
practical viewpoint. Structural vibrations and failures are 
often attributed to the periodic forces set up on members 
exposed to wind or other flow (smokestacks, pipe lines, struc- 
tural columns, to mention a few). In many cases it might 
be possible to destabilize the vortex shedding by addition 
of simple interference elements or by incorporating them in 
the original designs. In the case where one member is 
buffeted by the wake of another the same principle might 
j be applied. 

SPREAD OF VORTEX STREET 

It has been observed by most investigators that the spac- 
ing ratio A/Z is KArmdn’s value (0.28) close to the cylinder 
but increases rapidly downstream. The increase of hjl is 
mainly due to the increase of A, since I changes very little 
(fig. 21). In the stable range this is the result of viscous 
diffusion of the real vortices. 

Hooker (ref. 28) has made an interesting analysis. First, 
a real vortex has a core of finite radius; its center is the point 
of zero velocity and maximum vorticity. Hooker shows 
that in a vortex street, where the velocity field of the otlier 
vortices must be taken into account, the points of zero 
velocity and maximum vorticity do not coincide. The 
point of maximum vorticity is unchanged, but the point of 
zero velocity is farther away from the center of the street. 
As the vortex decays, the point of zero velocity moves 
farther out, its distance from the center of the street increas- 
ing almost linearly with time. Thus the spacing based on 
vorticity centers remains constant, while the spacing based 
on velocity centers increases linearly. Hooker's calcula- 
tion of the linear spread checks fairly well with some pic- 
tures taken by Richards (ref. 29) in the wake of an elliptical 
cylinder having a fineness ratio of 6: 1 and the major diameter 
parallel to the free-stream velocity. 
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However, the spread of the wake is not always observed 
to be linear. Among the different investigators there is a 
large variation of results, apparently dependent on the 
experimental arrangement. In Richards’ experiment the 
cylinder was towed in a water tank and the vortex patterns 
were observed on the free surface. 

In Kovasznay’s experiment the cylinder was mounted in 
a wind tunnel, the arrangement being similar to the one 
used here (see the section “Experimental Data”). On 
his plot of the streamlines at 7?=53 the downstream spread 
of the vortex street is parabolic rather than linear. It is 
possible to fit his results by a somewhat different applica- 
tion of Hooker’s idea, using decaying vortex filaments. 

Each vortex in the street is considered to behave like a 
single vortex filament carried along by the fluid, its decay 
or diffusion being the same as if it were at rest. The decay 
of such a vortex is described by a heat equation, whose 
solution is (ref. 30, p. 592): 

(BlO) 

wliere ^ is the tangential velocity at the distance r from the 
center and at the time t. The circulation is F. This is 
essentially a vortex with a ''solid'’ core and potential outer 
flow joined by a transition region in which the velocity has a 
maximum value. This maximum velocity is 

^*=0.72(r/27rr*) (Bll) 

and occurs at the radial distance 

r*=2.24(r0'/' (Bl2) 

Here r* is defined as the vortex radius. 

Thus the radius increases as and the maximum velocity 
decreases as r vortex street, the time t is replaced 

by the downstream distance x. Since the vortices move 
with the velocity rather than the dimensionless time 

Uo X . 

^ appropriate, where Z7* also varies downstream 

(see the section "Downstream Vortex Spacing"). 

Wlien a pattern of such vortices is superimposed on a uni- 
form flow, it is possible to calculate the velocity fluctuation 
at a point due to the pattern passing over it. 

Now the following hypothesis is added. It is assiuned 
that the vortex radius r* is equal to the width h of the street. 
Then the width of the street increases as 

A second result follows. The maximum velocity fluctua- 
tion (observed by a hot-wire, say) will occur on the line of 
vortex centers and will have the amplitude 

<1* (B13) 

that is, the hot-wire encounters instantaneous velocities 
varying from Z7* (because of vortex centers passing over it) 
to because of the fields of vortices on the other side 
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of the street). Relations (Bll) and (B12) then give the 
downstream behavior of the maximum fluctuation amplitude. 

The results may be summarized as follows: 

(a) Wake width A 

(b) The maximum amplitude of fluctuation u* occurs on 
the line of vortex centers (so there are two maximum points 
across the wake). 

(c) 

(d) «*=0.36 {Vl2wh). 

A comparison of the above predictions was made with 
calculations based on Kovasznay’s measurements which in- 
clude profiles of velocity fluctuation amplitude (cf. fig. 11) 
as well as the streamline plot. The following comparisons 
were obtained, item by item: 

(a') The time variation of hjd, determined from the vortex 
centers on the streamline plot is shown in figure 21. The 
parabola A/d=O.59(0— 6)‘'^ is shown for comparison. 

(b') The line of maximum velocity fluctuation lies slightly 
inside the line of vortex centers and is fitted by h*ld=0.53 

(e-6)'^l 

(c') The time variation of u* is also plotted in figure 21. 
(Actually Kovasznay’s maximum root-mean-square values 
u'„ are plotted, but these should differ from u* only by a 
constant factor.) The curve 'a'm/f^o=0.26(^ 6) is shown 

for comparison. The points could be fitted better, but the 
curve was chosen again to have the origin 6=6. 

(d') A comparison with (d) may be made by estimating 
the strength T of the vortices. Such a consideration, in fact, 
led to the present model, for it was found that the magnitude 
of the observed velocity fluctuations could be accounted for 
only by assuming that the radius of the vortex core is about 
equal to the width of the street. This observation had 
already been made by Fage and Johansen (ref. 8), for 
R«2X10^ If the free vortex layer is represented by a 
velocity discontinuity U=Uo to U=0, then the circulation 
is Uo per unit length and “the circulation” flows with the 
velocity (7„/2. On the other hand, the rate at which circu- 
lation enters one side of the street is mT, where T is the cir- 
culation per vortex.® Therefore 

_ Uo‘_ Uod 

2«i 2S 

For Kovasznay’s example, (SssO.lS, so T^iUcd. Then, 
comparing with (d), the maximum fluctuation in the initial 
part of the wake is 

^« 0 . 2 ^= 0.2 

Uo n, 

assuming h^d &t this low Reynolds number. The largest 
value of u'JUo in Kovasznay’s example is 0.14 at x/d=7, 
corresponding to u*IUo « 0.2. The order of magnitude of this 
estimate is quite sensitive to the size of the core relative to 

• The velocity at the outer edge of the layer is actually about 1.5 C7*, but experiments 
indicate that only about half the vorticity of the shear layer goes into individual vortices. 
Therefore, the value Uod/2S is a fair estimate. 


the width of the street; if the core is assumed to be much 
smaller, the calculated velocity fluctuations are much larger 
than those observed. Also, if the cores were very small 
compared with the width of the wake, four peaks instead of 
two would be observed in the profile of the velocity fluctua- 
tion amplitude. 
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TABLE I.— PROBABILITY FUNCTIONS 


u(t) 

p(f) 

F({) 

Nk 

c 

a 

Random function (Gaussian) 
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TABLE II 

RING DIMExNTSIONS 


Ring 

d, cm 

D 

d. 

Did 

1 

0. 168 

1. 59 

0.030 

9.5 

2 

.081 

.81 

.018 

10.0 

3 

.079 

.40 

.018 

5.1 


TABLE III 

VALUES OF STROUHAL NUMBER FOR VARIOUS TEST 
REYNOLDS NUMBERS 



R 

S 

Rd 

So 

89 

0.051 

450 

0.26 

96 

.052 

490 

.265 

103 

.052 

525 

.265 

128 

.057 

650 

.29 

153 

.060 

780 

.31 

182 

. 147 



215 

. 189 



302 

.204 



366 

.211 



455 

.212 
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Positive directions of axes and angles (forces and moments) are shown by arrows 


Axis 

Force 
(parallel 
to axis) 
symbol 

Moment about axis 

Angle 

Velocities 

Designation 

Sym- 

bol 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) 

Angular 

Longitudinal 

X 

X 

Rolling 

L 

Y >Z 

Roll 

fh 

u 

■n 

Lateral— 

V 

V 

PitcJiing 

M 

Z >X 

Pitr.h 

0 

iJ 

ft 

Normal 

Z 

z 

Yawing 

N 

X >Y 

Yaw 


X/ 

V) 

r 


Absolute coefficients of moment 

n 

”'~(LcS 
(pitching) 


(7— ^ 

^‘~qbS 

(rolling) 


^”~qbS 

(yawing) 


Angle of set of control surface (relative to neutral 
position), d. (Indicate surface by proper subscript.) 


4. PROPELLER SYMBOLS 


D 

P 

pID 

V' 

V. 

T 

Q 


Diameter 
Geometric pitch 
Pitch ratio 
Inflow velocity 
Slipstream velocity 


Power, absolute coefficient <7p= 


pnm^ 



a 

Speed-power coefficient 

Pn* 

rp 

V 

Efficiency 



n 

Revolutions per second, rps 

/ V \ 

Co——^— 


EfiPective helix angle =tan“' 

(27rm) 





5. NUMERICAL RELATIONS 


1 hp=76.04 kg-m/s=550 ft-lb/sec 
1 metric horsepower=0.9863 hp 
1 mph= 0.4470 mps 
1 mps=2.2369 mph 


1 lb =0.4536 kg 
1 kg=2.2046 lb 
1 mi= 1,609.35 m=5,280 ft 
1 m=3.2808 ft 




